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Abstract 

The notion of multidimensional quadrilateral lattice is introduced. It 
is shown that such a lattice is characterized by a system of integrable 
discrete nonlinear equations. Different useful formulations of the sys- 
tem are given. The geometric construction of the lattice is also dis- 
cussed and, in particular, it is clarified the number of initial-boundary 
data which define the lattice uniquely. 
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1 Introduction 



During the last three decades a new branch of mathematical physics, the 
theory of integrable systems, has been developed Besides its mathe- 

matical beauty, the theory has found also an enormous variety of applications 
in physics. 

A lot of attention has been dedicated also to discrete (difference) inte- 
grable systems [|ll01@|-[0]0-|l6|. They share with their differential coun- 



terparts many common properties, like for example an integrability scheme, 
constants of motion, soliton solutions, etc. 

To construct an integrable discretisation of the given soliton equation one 
can follow several different approaches. The standard ones are: 

i) a discrete version of the Lax pair (see, for example IQ, 0) 

ii) the Hirota method via a bilinear form [|10| [|lT| 



iii) extensions of the Zakharov - Shabat dressing method |T2 



iv) direct linearization using linear integral equations p!4|-||16||. 

Recently, another discretisation approach has been applied to those soli- 
ton equations which describe geometrically meaningful objects, like curves 
and surfaces. It is well known that many integrable PDE's allow for such an 
interpretation (see, for example [|T^-[2D| and references cited there). Given a 
geometric integrable system, the geometry provides its integrability scheme 
(the associated linear problem) and associates with such an integrable system 
other integrable ones, which describe various objects "living" on the same 
submanifold Jl^ (see, for instance, the connection between the integrable 
vortex motion and the Heisenberg spin model). 

It may be interesting to note that many of the important integrable PDE's 



were studied by the prominent geometers of the XlXth century ETJ |23]. They 



have found also tools which allow to construct explicit formulas for the im- 
mersed submanifolds, like the Darboux transformation pT[ |. 

In the process of discretisation of such systems, as a working principle, the 
discrete analogues of the relevant geometric properties must be found. In this 
way the discrete curves ||23| were " incorporated" into the soliton theory, and 
some classes of surfaces, including the pseudospherical, minimal, isothermal 
and constant mean curvature surfaces were discretised |2^-|26[. Recently a 
new example of integrable discrete surface, which contains in principle all 
the ones known before as a special reductions, has been found and related to 
the fully discrete Toda system (Hirota equation) 
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In this paper we discretise the geometrically meaningful integrable mul- 
tidimensional system which describes submanifolds parametrized by multi- 
conjugate coordinates. The associated differential equations were derived by 
Darboux |^ and rediscovered (and solved), in the matrix case, about ten 
years ago by Zakharov and Manakov p8|, using a d approach; for this reason, 
they are now called Darboux-Zakharov-Manakov (DZM) equations. More 
precisely, we do the following. 

i) We show that the proper geometric discretisation of the notion of (multi- 
dimensional) conjugate net is the notion of (multidimensional) quadrilateral 
lattice (MQL), meaning that all the elementary quadrilaterals of the lattice 
are planar. 

ii) We characterize MQL's in terms of a system of nonlinear integrable dis- 
crete equations in dimensions, > 3 (equations (0) of Section 2). 

iii) We show that the lattice is defined uniquely by ^^'^^^^ "initial" surfaces 
or, equivalently, by "initial" curves plus A^(A^ — 1) arbitrary "initial- 
boundary" data, which are functions of two variables. 

We remark that the equations which characterize the MQL have been 
recently obtained by Bogdanov and Konopelchenko [jl3| as a natural discrete 
integrable analogue of the DZM equations, but their geometric meaning was 
unknown. 

The MQL introduced in this paper allows one to generate, through a 
systematic geometrically meaningful reduction procedure, not only all the 
known geometric integrable lattices previously mentioned in this introduc- 
tion, but also new ones. For this reason, a systematic investigation of the 
reductions of the MQL has recently started |^ and the list of integrable 
geometric lattices already obtained through it includes the following. 

i) The fully discrete Toda lattice (Hirota equation) which, as it was shown in 
[p7[] , describes the Laplace sequence of 2-dimensional quadrilateral lattices. 

ii) The discrete analogue of the Lame equations for orthogonal nets, cor- 
responding to the case when the elementary quadrilaters are inscribed in 
circles. 

iii) Some symmetric reductions of the MQL including, for example, a discrete 
analogue of the (—1) element of the Nizhnik-Veselov-Novikov hierarchy of 
2 + 1 dimensional integrable systems. 

The integrability of the equations associated with the MQL was estab- 
lished in [O] using the d approach. A more geometric way of constructing 
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solutions of the MQL equations has been recently obtained in ||3T[, using a 
discrete vectorial variant of the Darboux transformations. These transforma- 
tions allow to construct, from a given MQL, a new one, which is, in general, 
topologically more complicated. In particular, dromionic and rational MQL's 
have been obtained. 

The paper is organized as follows. In Section 2 we give the definition of the 
multidimensional quadrilateral lattice and derive the underlying integrable 
system of discrete equations. In Section 3 we present different forms of the 
system. In Section 4 we discuss the initial value problem for the MQL's. 

2 Multidimensional Quadrilateral Lattices 

As it was mentioned in the introduction, the notion of quadrilateral lattice 
is the discrete generalization of the notion of conjugate nets. 

In the continuous case, conjugate coordinates on a surface are defined as 



follows [32|. Given a point on a surface and a second point which belongs 
to the coordinate line passing through the first one, the tangent planes to 
the surface at those two points intersect along one straight line. In the limit 
when these two points coincide, the direction of the straight line is the tangent 
direction to the second coordinate line (see Fig. 1). 



Figure 1. 

On a discrete level, a discrete surface can be defined as a mapping from 
1? to the M dimensional space M*^ 

. (1) 
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If we consider the two neighbouring "tangent planes" along the first coordi- 
nate line, defined by the triples {x,Tix,T2x) , {Tix,T^x,TiT2x) , where Tj is 
the shift operator in the i-th direction of the lattice, their intersection line 
coincides with the direction of the second coordinate line if and only if the 
elementary quadrilateral {x,Tix,T2X,TiT2x} is planar (see Fig. 2). 




Figure 2. 

Therefore the natural discrete analogue of the notion of 2-dimensional con- 



jugate net is given by the notion of 2-dimensional quadrilateral lattice [33 



Definition 1 A 2-dimensional quadrilateral lattice is a mapping x : — »• 
M^^ , M > 2 such that all the elementary quadrilaterals {x, Tix, T2X, T1T2X} 
are planar. 

A multidimensional conjugate net is characterized by the property that 
any surface made by varying two parameters forms a 2-dimensional conju- 
gate net with respect to these two coordinates. The discrete analogue of the 
above notion is therefore given by the notion of multidimensional quadrilat- 
eral lattice. 

Definition 2 By an A^-dimensional quadrilateral lattice we mean a mapping 
from an A^-dimensional integer lattice to the M-dimensional linear space 



X : Z 



N 



M> N 



(2) 
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such that all the elementary quadrilaterals with vertices 



{x,TiX,TjX,TiTjx} , i,j = 1...N , ij^ j , (3) 

are planar. 

Since the planarity of four points is preserved during the projective trans- 
formations of the ambient space, quadrilateral lattices are objects which ac- 
tually should be described whithin the projective geometry approach. In the 
above definition, therefore, it is very convenient to replace the linear space 
M*''^ by the projective space P^"'^: 

X -.Z^ , M > N . (4) 

To perform the calculations it is natural then to use homogeneous coordinates 
X G R^'^"*"^ \ {0} such that x = [x] G P''^ are the corresponding directions. 

In terms of the homogeneous coordinates the /m ear relations between the 
vectors representing vertices of the elementary quadrilaterals can be written 
system of ^'■^ discrete Laplace equations 

A A-x = iTiA,)D^^ + {T,Aj,)Dj^ + Cij^ , i^j, i,j = l,...,N , (5) 

where Dj = Tj — 1 is the discrete partial derivative, Aij and Cij = Cji are 
scalar functions of the lattice variable n = {ni, . . . , n^) G Z^: 

:Z^^M , ly^j, i,j = l,...,N . (6) 

The projective picture does not change if we multiply the homogeneous co- 
ordinates by a scalar function 

X i = ^x , p : ^ M \ {0} . (7) 

In fact, after such a gauge transformation, the new coordinates satisfy the 
Laplace equation with the coefficients 

Aj = 7^{Aijp~Djp) , z^j, i,j = l,...,N , (8) 
ijP 

Cij = -^{-D,Djp+{TiAi,)Dip+{T,Aji)Djp+Ci,p) . 
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Taking as gauge function (say) the last coordinate of x, i.e. p = x*^+^ (the 
non- homogeneous gauge), we "place" our net in the M- dimensional affine 
subspace of M^^+^ characterized by x*^"*"^ = 1. Then we get the special 
Laplace equations 

ADjx = (Ti^,,) Ax + {TjAj,)Dj^ , t^j, i,j = l,...,N , (9) 

where we can identify also the affine subspace with M^^, i.e. in the above 
formula, x = (x\ x*^) G M^^. 

For = 2 we have only one linear equation and no restrictions on its 
coefficients. The solution of the Laplace equation represents then a two- 
dimensional quadrilateral lattice (or discrete conjugate net) . 



Starting from iV = 3, we must repeate the two-dimensional construction 
in every pair of directions and, to make it possible, we should take into 
account the following compatibility condition between the Laplace equations 



DkAij — AijTjAjk + AikTkAkj — Ai^TkAij + Cij , (10) 
DkCij = DiCkj , i ^k^i . 

In the non-homogeneous gauge [Cij = 0) the above equations reduce to 
N{N -1){N -2) equations for the N{N - 1) coefficients Ay 

DkAij = AijTjAjk + A.ikTj.Aj.j — An^T^Aij , i ^ j ^ k i . (11) 

Therefore the system of equations (|1T]) is overdetermined for > 3 and 
determined only for = 3. 

The continuous limit of the equations and (|TT]) gives rise to the clas- 
sical Laplace equations 

and to the DZM equations for conjugate multidimensional nets PT 



= aijttjk + OLikCikj - aikttij , i ^ 2 ^k ^ i , (13) 

where /,^ = ^. 

We conclude this Section by pointing out that multidimensional quadri- 
lateral lattices are characterized by the set of integrable equations (O), and 
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can be therefore called integrable lattices. Their integrability scheme is gov- 
erned by the linear (discrete) Laplace equations (or (|^)) which express 
the planarity of the elementary quadrilaterals. 

In the rest of the paper we will refer to equations ( pT]) (or (p!OD) as to the 
multidimensional quadrilateral lattice equations. 

3 Alternative Formulations of the Multidi- 
mensional Quadrilateral Lattice Equations 

In this Section we discuss, for completness, other forms of the MQL equa- 
tions which will be usefuU in deriving integrable reductions [BO] and in the 



construction of the Darboux transformations |^ for the lattice. 

Equations (Fo|b), together with the symmetry condition Cij = Cji, imply 
the existence of a function C : M such that 

a, = D,D,C 1^3- (14) 



On the other hand from equations ([T(]|a) one can derive the following con- 

/ijj -|- i /ijfc -|- i 

which is automatically satisfied expressing Aij in terms of the "logarithmic 
potentials" 7^ : ^ R 

A., = ^ , . (16) 

Finally, the nonlinear system ([l0|) can be rewritten as a system of ^^^-^y^^-'^) 
equations 

(17) 

We remark that the last term vanishes when working in the non-homogeneous 
coordinates. 

To see the geometric meaning of the function 7, defined in ([T6|) let us 
note that, when taken as a gauge function, it removes all the corresponding 
coefficients Aij , {j = 1, . . . , N , j ^ i) from the Laplace equations (|^). 
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Of course, the definition of the functions 7^ depends on the particular 
gauge we start with. /^From now on, we fix the non-homogeneous gauge (|^) 
and define the functions 7^ with respect to that gauge. However, we still have 
the freedom of multiplying 7j by an arbitrary function of the single variable 

Ui. 

Let us rewrite the special Laplace equations (|^) using the functions 7^ 

DiDj^ = (t^^^ Ax + (^.^) ' ■ (18) 



Following Darboux [^, one can define suitably normalized tangent vectors 

X, = ^ (19) 

which possess the property that their variation in the second direction j ^ i 
is proportional to Xj only: 

D,Xi = (T,/3,,)X, , ly^j (20) 

where 

/3,. = ^ , . (21) 

The compatibility condition written in terms of the new fields f3ij takes the 
following form 

Dj/3,k = (3ij{Tj(3,k) , i^j^k^i. (22) 

As it was mentioned in the introduction, the matrix version of equations (|T7p- 
(|22D, for = 3, has been derived in |T^, through a d dressing approach, 
without a geometric interpretation, as an integrable discrete generalization 
of the DZM equations. 



4 Geometric Construction of the Lattice 

We conclude the Letter with some remarks about the constructability of the 
quadrilateral lattices. We discuss also the meaning of some objects intro- 
duced in the last Section. 
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Given two discrete curves xi(ni) and X2(n2) in M*''^ meeting in a point 
xi(0) = X2(0), and given two functions ^412 and ^421 defined on Z^, these data 
allow to construct uniquely the discrete quadrilateral surface. Equivalently, 
one can use instead of A12 and A21 the functions 71 and 72. Once the functions 
■ji, i = 1,2 are given, the vectors Xj are also defined in the points of the 
discrete quadrilateral surface according to formula (p!9|). 

The construction of the vectors Xj for a given 2-dimensional quadrilateral 
lattice (2DQL) is shown in the Fig. 3. For example, -D2X1 (the change of Xi 
in the second direction) is proportional to X2 only (the dashed line connecting 
the arrows of Xi and T2X1 is parallel to -D2X) . 



Let us consider now three discrete curves Xi(ni), X2(n2) and X3(?7,3) in 
(M > 3) meeting in a point xi(0) = X2(0) = X3(0). For any pair of 
curves: Xj(nj), Xj(nj), assigning two arbitrary functions Af^^ : — M we 
construct a 2-dimensional (initial) QL. It turns out that, starting from these 
three initial discrete surfaces, one constructs uniquely the three dimensional 
quadrilateral lattice. The construction is based on the simple fact that three 
different, non parallel, planes in a 3-dimensional space intersect, in general, 
in a single point. It may happen that the intersection point is at infninity. 




Figure 3. 
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but we still have non ambiguous construction looking at the lattice in the 
projective space P*^ = M*^ U P*^-^ 

In particular, the point T1T2T3X is the point of intersection of the three 
planes (Tix, TiTsx, TiTsx), (Tgx, TiTax, T2T3X) and (T3X, T1T3X, TsTgx) in 
the three dimensional subspace (x, Tix, r2X, Tax) of M-^ (see Fig. 4). 




Figure 4. 

The above construction involves only linear operations and can be pro- 
longed to build, out of the three initial quadrilateral surfaces, the whole 
3-dimensional lattice. 

Fixing the vectors Xj (or the functions 7^) on the initial curves allows to 
find their values in the points of the initial discrete quadrilateral surfaces, 
and then in the points of the whole lattice. Because of the planarity condition 
there is no contradiction in the construction of the vectors Xj. For example, 
the vector T2T3X1 is obtained by the intersection of the line (T2T3X, TiT2T-^x.) 
with the plane parallel to (x, T2X, T3X) passing through the arrow of Xi (see 
Fig. 5). 
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Figure 5. 

To summarise, in order to construct uniquely the 3-dimensional quadrilateral 
lattice, one has to give three arbitrary intersecting quadrilateral surfaces or, 
equivalently, three arbitrary intersecting curves plus six arbitrary functions 
of two discrete variables 

A^\ni,nj) ,Af^{m,nj) , l<i<j<3 . (23) 

The construction outlined above allows then to build uniquely the whole 
3-dimensional lattice or, equivalently, the functions 

Aij{ni,n2,n3) , i j , i,j = 1,2,3 (24) 

which solve the MQL equation and satisfy the following boundary conditions: 

Ai2(ni, 712, 0) = 74^2^(^1,712) , ... ,^32(0, 712,713) = ^^2^ (712,713) • (25) 

If, instead of the Aj^'s, we use the data 7j's, then our arbitrary initial func- 
tions are 7^^(711,712) , 7i°Vi>'^3), 72°Vi>'^2), 72°V2,»^3), Ts^Vij^-s) and 

Let us apply the same procedure to the case = 4. Given four discrete 
curves and six, constructed from them, arbitrary initial 2DQL's, they give rise 
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to the four unique 3DQL's x(ni, n2, n^, 0), x(ni, n2, 0, n^), x(ni, 0, ris, n^^) and 
x(0, 77-2, ns, in the way described above. Then the point x(l, 1, 1, 1) can 
be constructed, for example, out of the surfaces x(l, n2, na, 0), x(l, n2, 0, 714) 
and x(l, 0, 71,3, n^), i.e. it is the intersection point of the three planes T1T2F34, 
T1T3P24 and T1T4P23 in the three dimensional subspace T1P234 of M^. Here 
Pij {i 7^ j) denotes the plane passing through the points x, TjX and T^x, 
i.e. Pij = Pij{y^) = (x,rjX,rjx), e. g., T1T2P34 = {TiT2X,TiT2T-iyi,TiT2T4x); 
analogously Pijk — (x, T^x, TjX, Tfcx) (all the indices are distinct). The same 
point could be constructed, however, also in three other ways. In the follow- 
ing we show that the results coincide. 

First, we note that the plane T1T3P24 is the intersection of the two 
three dimensional subspaces T1P234 and T3P124 in the four dimensional space 
^1234 = (x, TiX, r2X, Tax, r4x). Similarly riT2P34 = T1P234 n r2Pi34 and 

T1T4P2S — T1P234 n J4P123. 



T2T3T4X 



T3T4X 



T1T3P24 



-A 




Tix 



Figure 6. 



The point TiT2T3T4:k is then the unique intersection point of the four three 
dimensional subspaces TiP^i^^, i = 1,...,4 (the symbol l.iA denotes the 
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sequence of the natural numbers from 1 to 4 with the ith element removed) 
of the four dimensional space P1234; consequently, it is the intersection point 
of the four triplets T^TiP^ -- k = 1, . . . ,4 (where the index i, i = 1, . . . ,4, 
i ^ k, indicates each element of the A;th triplet) of the two dimensional spaces 
in the corresponding three dimensional spaces T^P^j^^ 

4 4 

Pu3A^TiT2nn^ = f]T,p,u= fi nr.p.u^ . (26) 

i=l i=l,i^fc 

The same argument can be used to prove the compatibility of the con- 
struction for an arbitrary dimension N of the lattice. In the natural notation 
inherited from the example = 4 

N N 

Pi...A.9Ti...TA.x=f|r,Pu.iv= n TkT.Pu.k.N , k = l,...,N 

(27) 

which allows to apply the mathematical induction. 

To summarise all the geometric considerations of this Section, we conclude 
it stating the general initial boundary- value problem for MQL's. 

In order to construct uniquely the N -dimensional quadrilateral lattice, one 
has to give arbitrary intersecting quadrilateral surfaces or, equiv- 

alently, N arbitrary intersecting curves (in general position) in plus 
N{N — 1) arbitrary functions of two discrete variables 

A^\n„n^) , Af^{ni,nj) , 1 < z < j < iV . (28) 

The linear construction outlined above allows then to build uniquely the 
whole N -dimensional lattice or, equivalently, the functions 

Aj(ni,n2 . . . ,n7v) , i^i , i,j = l,...,A^ (29) 



which solve the MQL equation (1J_) and satisfy the following boundary con- 
ditions: 

Aj(0, . . . ,ni, . . . . . . ,0) = AfMni.Uj) , (30) 



Aji(0, . . . ,ni, . . . . . . ,0) = Af^{ni,nj) , l<i<j<N . 
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Note that the number of the arbitrary functions of two variables enter- 
ing into the general solution of the MQL equation agrees with that of the 
continuous case lET . 



Remark All the geometric considerations of this Section can be reformulated 
in terms of linear equations (for example, a 3-dimensional subspace of a 5- 
dimensional space is equivalent to a system of (5-3) linear equations for 5 
unknowns) and the corresponding theorems about their solutions. 

5 Final Remarks 

It is important to emphasize that the MQL is, by construction, an integrable 
lattice. Indeed, the linear constraints (||) (the planarity constraints) pro- 
vide a well defined way to construct the lattice; therefore they charaterize 
completely the MQL and its integrability properties. 

We are convinced that the simple geometric meaning of the MQL equa- 
tions presented in this paper will allow to find for them many applications in 
all the branches of physics in which studies of multidimensional lattices are 
of importance. 

We would like to point out that the projective picture used in this paper 
should not be considered as a mathematical curiosity. It allows to treat the 
points of the MQL which belong to the hyperplane at infinity on an equal 
footing with the "usual" points of the affine space. The presence of the points 
at infinity corresponds to singularities in the solution of the MQL equation 
(pUf). It may be interesting to develop this point of view in connection with 
the singularity confinement property of integrable discrete systems ||34|| . 
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